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Abstract We discuss an autonomous motor based on a Brownian particle 
driven from thermal equilibrium by periodic in time variation of the internal 
potential through which the particle interacts with molecules of the sur¬ 
rounding thermal bath. We demonstrate for such a motor the absence of a 
linear response regime: The average driving force and drift velocity are shown 
to be quadratic in both the frequency and amplitude of the variation. The 
adiabatic approximation (of an infinitely slow variation) and the leading cor¬ 
rection to it (linear in the variation’s frequency) both lead to zero drift and 
are insufficient to describe the motor’s operation. 
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1 Introduction 

Consider a Brownian particle moving in one dimension and interacting with 
molecules from the left and right through microscopic potentials Ui and Ur 
which are of different shapes and/or ranges. Contrary to uncultivated intu¬ 
ition (but in agreement with thermodynamics), in thermal equilibrium such 
intrinsic microscopic asymmetry does not cause a net drift of the particle: 
Although a molecule, say, from the right interacts with the particle via a 
stronger force than a symmetrically positioned molecule from the left, the 
average forces exerted on the two sides of the particle, calculated with the 
equilibrium Boltzmann distribution, have exactly the same magnitudes and 
completely compensate each other. On a deeper level, this cancellation is en¬ 
sured by detailed balance symmetry [T] which may (or may not) be broken 
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when the system is out of equilibrium. Now suppose that the particle is active 
in the sense that it is equipped with an internal mechanism which modulates 
the static potentials Ui and Ur in a periodic-in-time manner 

Ur,it)=Ut)Ur., ( 1 ) 

e.g. with harmonic modulation functions 

Cc((^) = 1 + «c( sinwat (2) 

with amplitudes 0 < Uq. < 1. Here and below the subscript a = {l,r} refers 
to the left and right sides of the particle. Such modulation may result from 
periodic conformational changes of internal degrees of freedom, which are 
very common for variety of biological macromolecules like proteins, ribo¬ 
somes, and viruses. Now when the particle is out of equilibrium, one may 
reasonably expect that the asymmetry of static potentials Ua, and/or the 
modulation amplitude Oq, and frequencies Wq, may result in the particle ac¬ 
quiring a nonzero average velocity and serving as a Brownian motor. As for 
many other machines operating under non-equilibrium conditions, the ex¬ 
plicit evaluation of the drift force and velocity is a non-trivial problem, and 
even the direction of the drift may be not easy to guess. There are also a 
number of peculiar aspects of this model that we believe make it worthy to 
discuss. 

Unlike many mesoscopic machines driven from equilibrium by external 
means (e.g., due to a contact with thermal baths of different temperatures, 
periodically in space and time modulated temperature, external potential, 
light, etc.), our motor is autonomous and resists the thermalization by means 
of the internal mechanism. Among other types of autonomous motors stud¬ 
ied in recent years are granular Brownian systems [2], chemically powered 
motors driven by asymmetric catalytic activity [5], and Brownian informa¬ 
tion machines [4]. An autonomous active motor with an internal anchoring 
mechanism was studied in [S]. 

Driven by the oscillation of a microscopic parameter(s), our motor cannot 
be described within empirical approaches, e.g. those based on the standard 
Langevin or Fokker-Planck equations with time dependent external param¬ 
eters [H[7j or modified with the energy depot (“negative friction”) terms [5] . 
Only a few microscopic models of rectified Brownian motion are presented so 
far in literature, and most of them concern the regime beyond the weak cou¬ 
pling to the thermal bath. In such cases, the Fokker-Planck [annuiiiiiiiiiii 
[B] and Langevin [BUIS] equations involve additional terms (of higher orders 
in a weak coupling parameter), and must be supplemented with additional 
fluctuation-dissipation relations which cannot be established phenomenolog¬ 
ically. In contrast, for the present model a systematic driving force will be 
shown to emerge already in the leading order in the weak coupling param¬ 
eter. Still, an explicit expression of the driving force cannot be constructed 
empirically and requires a microscopic evaluation. 

One interesting property of the model is the absence of the linear re¬ 
sponse regime with respect to the modulation frequencies ui considered as 
a perturbation parameter. For low w, the driving force and drift velocity 
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depend on frequency as Neither the adiabatic approximation (asymp¬ 
totically slow modulation), nor the leading perturbational correction to it 
are sufficient to account for directional motion of the particle. A similar be¬ 
havior was found previously for Brownian ratchets driven by temperature 
oscillations [BimiTT] . Other systems driven from equilibrium by variation of 
external parameters may show also qualitatively different scenarios ranging 
from directional transport driven by adiabatically slow variations [TOII^ to 
the absence of net transport beyond the adiabatic regime mm- Whether 
such diversity can also be observed in the family of autonomous machines 
like ours is yet to be explored. 

It might be relevant to note that the absence of the linear response regime 
was previously observed also for motors with non-varying parameters but 
driven by the coupling with two baths with non-equal temperatures mm 
m. In that case, the prediction that the drift velocity must be an even func¬ 
tion of a perturbation parameter (the temperature difference) can be envis¬ 
aged from simple symmetry arguments. No similar arguments are apparently 
available for the present model. 

Another remarkable feature of the model, which also has counterparts 
among motors driven by variation of external parameters [niiiH] , is a non¬ 
monotonic dependence and - for certain regimes - the sign reversal of the 
drift velocity as a function of oj (see Fig. 3 below). However, for our model 
such behavior is probably of only academic interest since it occurs when the 
period of modulation u}~^ is unrealistically short - of order or shorter than 
the collision time r. Since the latter is usually the shortest characteristic 
time, we shall restrict our theoretical discussion to the low frequency limit 
w T <C I. Numerical simulation will be used to confirm theoretical predictions 
and to extend the results beyond the low frequency domain. 


2 Model and simulation 

We consider a Brownian particle of mass M immersed in the thermal bath of 
temperature T comprised of ideal gas molecules of mass m <C M. Velocities of 
molecules before collisions with the particle are distributed with the Maxwell 
distribution 

where vt = \/kB T/m is a thermal (mean-squared) velocity of a molecule. 
The motion of the particle and molecules of the bath is assumed to be one¬ 
dimensional. 

Dynamics of the particle’s internal degrees of freedom, are governed by 
a certain internal built-in mechanism whose specific design is immaterial for 
our purposes. When this mechanism is turned off the particle is “passive”, i.e. 
behaves as a conventional Brownian particle interacting with molecules of the 
bath with time-independent forces. In order to facilitate analytic calculations, 
we assume that in the passive regime the bath molecules, unless their energy 
is too high, interact with the particle with constant repulsive forces of a finite 
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Fig. 1 The solid line is the potential energy for the interaction of the motor (large 
circle) in the passive regime and a molecule of the bath (small circles) as a function 
of position of the latter. When a molecule is inside the left (right) interaction zone 
X > Xi (x < Xr), it experiences a constant repulsive force of the amplitude fi (fr)- 
It is assumed that temperature is sufficiently low, so that only a negligible number 
of molecules with energy much higher than ksT experience the upper nonlinear 
part of the potential. Dashed lines represent a periodic in time modulation of the 
potential for the active regime. 


range, see Fig. 1. Namely, the potential energy of interaction of the particle 
with bath molecules on its left reads 

Ui= fiY,i^i-Xi)e{x,-Xi). (4) 

i 

Here Xi are coordinates of bath molecules, fi is a positive constant of the force 
dimension, d{x) is the step-function. The coordinate Xi is associated with 
the left side of the particle and determines the boundary of the interaction 
zone for molecules coming from the left: molecules outside the interaction 
zone, Xi < Xi, do not interact with the particle, while every molecule inside 
the zone, Xi > Xi, exerts on the particle the same force fi. Similarly, the 
potential energy of interaction with molecules on the right of the particles is 

Ur = - Xr) 0{Xr - X, ), (5) 

3 

where the positive constant fr is the amplitude of the constant repulsive force 
exerted on the particle by a bath molecule from the right when the latter is 
in the right-hand side interaction zone, Xj < Xr- 

To avoid complications related to the overlapping of left and right in¬ 
teraction zones, one can assume that closer to the particle’s core the linear 
potential is replaced by a sharper (diverging) one, see Fig. 1. This, however, 
is of no consequence as soon as the crossover from the linear to nonlinear 
potentials occurs at the energy much higher than ksT: only a negligible 
fraction of molecules with velocities v ^ vt would feel the nonlinear part of 
the potential. 

For this model the fluctuating Langevin force F{t) and its correlations can 
be readily evaluated analytically (see section 3 in [23]). To the lowest order 
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in the mass ratio parameter A = \Jm/M and for the time scale much longer 
than the collision time (defined by Eq. (14711 below) the Langevin equation 
for the particle’s velocity V has the standard form 

M^ = -^V + F{t). ( 6 ) 

A possible asymmetry of the microscopic forces, fi ^ fr, does not show up 
in this equation. Although the correlation function of the Langevin force 
{F{0)F{t)) does depends on /; and fr, this dependence disappears after the 
integration over time. As a result, the damping coefficient 7 does not depend 
on parameters of microscopic dynamics and takes the form 

I r°° [2 

7 =-— / {F{0)F{t)) dt = -nmvT, (7) 

kbT Jo V tt 

where n is the concentration of bath molecules. Also, if one writes E as a 
sum of forces on the left and right sides of the particle, F = Fi + Fr, one can 
show that 


{Fi) = -{Fr)=nkBT, ( 8 ) 

so the net Langevin force is zero-centered, {F{t)) = 0. In accord with thermo¬ 
dynamics, the microscopic asymmetry does not bias equilibrium Brownian 
motion. 

Our goal is to generalize the above microscopic model of passive Brownian 
motion for the case when the particle is “active”, that is capable to variate 
the potentials Ui and Ur according to Eq. O- Eor the active particle the 
potential energy of interaction with the bath is given by expressions similar 
to that for the passive regime 

Ui(t) = fi(t) - Xi) e{xi - Xi) 

i 

Ur{t) = -frit)J2ix, - Xr) 0(A, - Xj), (9) 

3 

in which the static force magnitudes fi and fr are now replaced by time- 
dependent ones. 


fa{t)=^ait)fa, ( 10 ) 

where ^ait) is given by ([ 2 ]), ^^(t) = 1 -|- sinwct (to minimize the number 
of parameters we assume no phase shift). Clearly, such a variation requires 
an external input of energy, which is not explicitly reflected in the model. 

The modulation frequencies Wq will be assumed to be small compared to 
the characteristic collision time, while the amplitudes Uq are not necessarily 
small and may take any values from the interval (0,1). The values Oa = 1, 
though not altogether meaningless (may correspond to a permeable particle, 
see [23]), are not considered. The variation of the particle-bath interaction 
supports the particle in a nonequilibrium state, which is the first condition of 
the drift. The second condition - the break of the spatio-temporal symmetry 
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time 

Fig. 2 Average dis plac ement {x{t)) and velocity {v{t)) of the motor as functions 
of time for regime GH) with parameters ujr = oJi = 0.1, ar = ai = 0.5, and 
S = fr/ fi = 3. Numerical experiment data (solid lines) are averaged over about 
2 • 10® tra ject ories. Theoretical curves (dashed lines) are solutions of the Langevin 
equation CT which corresponds to the average fluctuating force {F{t)) given by 
Eq. (l54l) . Units are defined by Eqs. dSSl). 


- can be arranged by assigning different sets of values {/a,act,WQ;} for the 
left and right sides of the particle. 

We study this model both theoretically and using numerical simulation. 
The latter is performed for the mass ratio = m/M = 0.01 with the 
standard molecular dynamics scheme with the only difference that instead of 
using periodic boundary conditions, we generate in the beginning of each sim¬ 
ulation run a very large thermal bath of noninteractive molecules (see [53] for 
details). The simulation shows that for unequal sets {/;, o/, w;} and {/r, a^, Wr} 
the particle develops non-stationary drift velocity, whose direction and value 
depend on fa,o,a,L>Ja in a rather subtle way. 

As a showcase example let us consider the case when asymmetry is due to 
non-equal magnitudes of the static (passive) forces, // ^ fr- Specifically, sup¬ 
pose that the particle, while in the passive regime, interacts with a molecules 
on the right via a stronger static force than with a molecule from the left (as 
in Fig. 1), while the modulation parameters for the left and right sides are 
the same, 

fl < fr, UJl=UJr= OJ, ai= Gr = Q. (11) 

For this case, when the modulation frequency is sufficiently low, the sim¬ 
ulation shows that the particle drifts to the right, that is in the direction 
of the steeper slope of the internal potential, see Fig. 2. Perhaps somewhat 
counter-intuitive, this behavior is in agreement with a theory which will be 
developed in sections to follow. Note that for the given values of parameters 
the drift is rather small and becomes visible only after taking average over 
many particle’s trajectories. A single trajectory looks like a random path of a 
passive Brownian particle and shows no visible bias. All experimental curves 
in this paper represent data averaged over about 10® trajectories. 
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Fig. 3 Drift velocity of the mo tor averaged over time as a function of the modula¬ 
tion frequency uj for regime (HB with parameters ar = ai = 0.5, and 5 = fr/fi = 3. 
Each experimental point corresponds to an ensemble-averaged trajectory similar to 
that in Fig. 2. 


Another interesting feature suggested by simulation data is that for low ui 
the drift velocity increases with frequency as This can be interpreted as 
the absence of the linear response regime with respect to w as a perturbation. 
As frequency is getting higher, the frequency dependence of the drift velocity 
becomes nonmonotonic, and the inversion of the drift direction occurs, see 
Fig. 3. 

According to m. molecules in the left and right interaction zones of the 
motor experience the forces of magnitudes 


fai^^ — fa “t“ fa ^a sinC^Q-t. 

(12) 

One observes that for the regime defined by Eq. dTlT) both static and time- 
dependent parts of these expressions are different for the left and right sides 

of the particle. It is natural to ask if the drift still occurs 
affects only the static parts. 

when asymmetry 

fl ^ fr fl^l — fr ^r: — ^r- 

(13) 

or only dynamic parts, e.g. 


fl ~ fr Gl (Xr^ Gll — l^r- 

(14) 


The simulation shows the drift in both these cases too, but with no sign 
reversal and in opposite directions, see Fig. 4. This suggests that the reversal 
of the drift direction for regime (ED can be interpreted as a result of the 
interplay or interference of two regimes described by Eqs. ED and ED- 
Another interesting regime of the motor’s operation is when the drift is 
induced by unequal modulation frequencies (one frequency may be zero), 
while the pair of other parameters is the same for both sides of the particle. 


fl — fr a/ — Gr: 


UJl ^ COr- 


(15) 




0.02 



frequency 

Fig. 4 Drift velocity of the motor averaged over time as a function of the modu¬ 
lation frequency u. The upper curve corresponds to the regime (I13II (asymmetric 
static part of inte raction) with fr/ fi = 3 and fi ai = fr ai — 0.5. The bottom curve 
is for regime (1141) (asymmetric dynamical part of interaction) with ar = 0.25 and 
Or = 0.5. 


In this case, if one or both frequencies are low (much shorter than the collision 
time), the particle drifts in the direction of the side with lower modulation 
frequency (two upper curves in Fig. 5). On the other hand, if both frequencies 
are high, the particle systematically moves in the direction of the side with 
a higher frequency (a bottom curve in Fig. 5). 

In the following sections we shall focus on a theoretical description of 
the model. We shall not try to cover the whole rich phenomenology which 
the model shows in simulation, but rather restrict ourselves to developing a 
perturbation approach for the low frequency domain. 


3 Theory: basic relations 

We assume that in the active regime the particle is still described by the 
Langevin equation with the damping coefficient approximately the same 
as for the passive regime, i.e. given by but with a non-zero centered 
fluctuating force, {F{t)) ^ 0. The goal is to evaluate microscopically {F{t)) 
and then, solving the averaged Langevin equation 

Mj^{Vit))=-^{V{t)) + {Fit)), (16) 

find the average velocity {V{t)) and trajectory {X{t)) of the motor. 

An obvious drawback of this approach is that it neglects the influence 
of the active part of the fluctuating force F{t) on the damping coefficient 
7 . Since the two quantities are related by a fluctuation-dissipation relation, 
this approximation cannot be entirely consistent. Yet it is clear that for suffi¬ 
ciently small uja or/and Oq an “active” correction to the dissipation constant 
is small compared to the value of the latter for the passive regime, and thus 
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time 

Fig . 5 Average displacement {x{t)) of the motor as a function of time for regime 
m when asymmetry is due to unequal modulation frequency, toi ^ lOr. The upper 
line is for uji = 2.5 and cur = 0, the middle line is for uji =0.1 and otj. = 0, and the 
bottom line is for uui = 1.1 and cJi- = 1. For all lines the modulations amplitudes are 
the same ai = ar = 0.5. The oscillatory character of the upper line is imperceptible 
on the figure’s scale. 


should produce a little effect. We shall see that the comparison of the theory 
with simulation results supports this intuition. 

As known from the microscopic theory of Brownian motion, the Langevin 
equation corresponds to the lowest order approximation in the mass ratio 
m/M, in which case the fluctuating force F{t) can be evaluated neglecting 
the particle’s motion. Setting F = Fi+Fr, let us consider the force Fi exerted 
on the particle, fixed in space, by molecules coming from the left. For the 
linear potential ®, each molecule in the left interaction zone x > exerts 
on the particle the same time dependent force Then the total force on 
the left side is given simply by the product 

Fi{t) = Mt)Ni{t), (17) 

where the Ni{t) is the number of molecules in the left interaction zone at a 
given time. 


Ni{t) = / dv dxf{x,v,t). (18) 

9-00 Jxi 

Here f{x, v, t) = — Xi)5{v — Vi) is the microscopic density of molecules 

in the position and velocity space. For simplicity we extended the right border 
of the left interaction zone to infinity. For low enough temperature the effect 
of such approximation is negligible. 

It is convenient to define the collision time Ti(v,t) as a time spent in the 
interaction zone x > Xi hy a. molecule which enters the zone with velocity 
V and leaves the zone at time t. The collision time Tr.{v,t) for molecules 
hitting the particle’s right side is defined in a similar way. The dependence of 
Ta {v, t) on time reflects the dynamical nature of the interaction potentials ® . 
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Fig. 6 Solid line circles: at time t molecule two is just before the collision entering 
the left interaction zone x > Xi, and molecules one is just after the collision leaving 
the interaction zone. Dashed line circles show same molecules at time t — ri. 


Of course, for the passive regime the collision times depends on molecule’s 
pre-collision velocity only. 

Consider two molecules with the same pre-collision velocity n > 0 at the 
border of the interaction zone Xl at time t. Molecule one is just after the 
collision, departing the interaction zone. 


xi{t)=Xi, wi(t)<0, (19) 

while molecule two is just before the collision, entering the interaction zone, 

X2{t) = Xi, V2{t) = v>Q, ( 20 ) 

see Fig. 6. According to the definition of Ta{v,t), at earlier time t — Ti(v,t) 
both molecules had velocity u > 0 and coordinates 

Xl{t-Tl)=Xl, X2{t-Tl)=Xl-VTl. (21) 

For the low modulation frequency it is reasonable to assume that molecules 
with the same initial velocity do not bypass each other in the interaction 
zone. Then it is clear that all molecules with a pre-collision velocity v > 0 
which are in the interaction zone at time t, at the time t — Ti were within the 
interval 


[Xi{t - Tl),X2{t - Ti)] = [Xi - VTl,Xi] (22) 

and moving to the right with velocity v, see Fig. 6. Thus for the number of 
molecules in the interaction zone at a given time, instead of (HHi) one can 
write 


poo pXi 

Ni{t) = / dv dx f{x,v,t - n). 

Jo j Xi—VTi 


(23) 


The advantage of this form is that it involves integration over coordinates 
and velocities of molecules before collisions with the particle, in which case 
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the microscopic density is simply the Maxwell distribution /m(w) multiplied 
by the average concentration of bath molecules n, 


{f{x,v,t)) =nfM{v). 


(24) 


Then from (I2dll for the average number of molecules in the interaction zone 
one obtains 


pOO 

{Ni{t))=n / dv fM{v)vTi{v,t), 
JQ 


(25) 


and the average force on the left side of the particle {Fi{t)) = fi{t) {Ni{t)) 
takes the form 


poo 

{Fi{t)) = n fi{t) / dv fM{v)vri{v,t). 

Jo 

Similarly, the average force on the right side of the particle is 

{Fr{t)) = -n fr{t) f dv fM{v)vTr{v,t). 

J —oo 


(26) 


(27) 


To advance further, we need to evaluate collision times Ta{v,t) for the 
active regime, which will be the focus of the next section. Meanwhile, it is 
instructive to consider how the general expressions (ESI) and EZl) work for the 
passive regime, when fa{t) = fa = const. In this case the collision times do 
not depend on time and we shall denote them as t^(v). Consider a molecule 
which enters, say, the left interaction zone at the moment t = 0 with the 
initial velocity v > 0. While the molecule is still in the zone its velocity 
evolves as 

v(t) = V — — t. (28) 

m 


The collision time corresponds to the moment when the molecule leaves 
the zone with velocity —v. The equation n(r°) = —v gives 




2mv 

~T 


(29) 


Similarly, the right side collision time is 


r°(v) = 


2mv 


(30) 


Substitution of these expression into Eas. (E5)) and (l?7l) gives the result (|H]) 
of the elementary kinetic theory 

{Fi)=-{Fr) = FQ = nm{v‘^) = nkBT. (31) 


As expected, for the passive regime the asymmetry fi ^ fr does not 
induce a net driving force on the particle. We can now interpret this no- 
go result as follows: Each molecule in the a-th interaction zone pushes the 
particle with a force of magnitude fa, but the average zone’s population 
(Na), according to (1251) . is linear in the collision time Tq, which in turn is 
inversely proportional to fa, Ta ^ fa- As a result, for the average force 
{Fa) fa i^a) the dependence on fa is canceled out. 
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4 Collision time for active regime 

Above we defined the collision time Ta{v,t) as functions of the pre-collision 
velocity v and after-collision exit time t, which corresponds to the moment 
when a molecule is leaving the interaction zone. With these functions one 
can express the forces {Fa{t)) on two particle’s sides in simple forms (1^ 
and (EZl). The disadvantage of functions Ta{v,t) is that they are difficult to 
evaluate directly. Let us define collision times t*{v, t) whose first argument is 
still velocity before the collision, while the time argument refers now to the 
moment when a molecule enters the interaction zone. Since both arguments 
of T*(u, t) refer to the same pre-collision moment, these functions are easier to 
evaluate. The two functions Ta{v,t) and T*{v,t) are related by the equation 

T{t) =T*{t-T{t)). (32) 

Here and for the bulk of this section we omit for brevity the argument v 
which is assumed to be the same for all quantities, and also suppress the 
left/right index a. We shall restore a in the final expression for Tq,. 

In order to find an explicit expression for r in terms of r* we approximate 
the right hand side of Eq. (l32ll by the first three terms of the Taylor expansion 
about t, 


(33) 

where primes denote time derivatives. Since (r*)' ^ w and {t*)" ~ the 
above relation is of second order in w. To the same order, the solution of Eq. 
(155)) has a form 

r = r* - (r*)' r* + 1 (r*)" (r*)^ + [(r*)f r*, (34) 

which produces a desirable explicit expression of r in terms of r*. Our goal 
now is first to evaluate r*, and then using (1551) to find r. Substitution of r 
into Eqs. (ESI) and (1271) will give us a perturbation expression for the force 
on the particle to second order in w. 

Consider a molecule which enters the left interaction zone x > Xi with 
velocity uq > 0 at the moment t = tg. Setting for a moment A; = 0, the 
equation of motion and initial conditions read 

mx”it) = -f{t) = 

x{to) = 0, x'(to) = Vo, (35) 

where ^(t) is given by (15|). The solution of the initial value problem (1551) is 
convenient to write as a function of time t elapsed since the moment to when 
the molecule enters the zone: 

x{t) = Uo t ^ -f 
2m 

^ ( sina;(to + “ sina;to — cosa;io) ■ (36) 

mcu^ V / 
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Recall that for the passive regime, the collision time r* can be determined 
from the equation v{t*) = —v, since the speed of a molecule before and 
after the collision is the same. That is, of course, not so for a time-dependent 
potential. For the active regime the collision time T*(tio,to) should be found 
as a nonzero solution of the equation 


a;(T*)=0. (37) 

We wish to find an approximate solution of this equation to order (an 
appropriate dimensionless small parameter is introduced by equation 
below). Using in (1361) the truncated expansion 


sina;(to + t) ~ sinwto + cosuito (wt) — - sinwto (wt)^ 

— - coswto + 777 sincuto 
6 24 

one obtains to order 

, , , , f a Lu cos ojtn o f o sin utn a 

6 m 24 m 

Here the first term is of zero order in u 


(38) 


(39) 


Xad{t) = Vot - (40) 

2 m 

and corresponds to the adiabatic approximation which completely neglects 
the change of the potential during the collision. As we shall see below, the 
adiabatic approximation x(t) « Xad(t) is not sufficient to account for the 
motor’s drift. The two last terms in (15^ are of first and second order in w, 
and take into account the dynamical nature of the potential. 

Substitution of dSH) into and resetting {vo,to) —>■ {v,t) gives for 
the collision time T*{v,t) a cubic algebraic equation which we write in the 
following dimensionless form: 

^ coswt ® (^) = 

where, recall, = t^{v) = 2mv/ f is the collision time for the passive regime. 
Solving this equation perturbatively to order one obtains (see Appendix): 

T* = T°|^“^(t) — ^ (uJT^) COSUjt 

+ ^ (wr°)^ sinujt + (wr°)^ cos^a;t|. (42) 

12 9 J 

Recall that T*(v,t) is a function of a time when a molecule enters the 
interaction zone, while the expressions (E51) and (l?fl) for the average forces 
involves the collision time t(u, t) as a function of the time when a molecule 
leaves the zone. To second order in to the relation between r* and r is given 
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by Eq. ([34]). Substituting (1421) into (IMl) . neglecting terms of order higher than 
and restoring the right/left index a we obtain 

T-a = T-° ^ (War°) (t) COSUJat 

+ SmUJat+ (War°)^ COS^W„t|. (43) 

To zeroth order in uj this expression gives 

Taiv,t)PiT°{v)^-\t) = (44) 

This is just the expression for the collision time for the passive regime, 
given Eqs. (P^)) and O, in which the static force magnitudes fa are replaced 
by dynamic ones fa{t) = fa^a{t)- Thus the expression (H41) corresponds to 
the adiabatic approximation and, as one can check, can be derived as a 
solution of the equation Xad{T*) = 0 with Xad{t) given by (HHll . 

Three last terms in the right hand side of Eq. (H51) describe corrections to 
the adiabatic approximation up to order These terms behave differently 
being time-averaged over the period of modulation: the term linear in w 
vanishes, whereas the terms quadratic in uj do not vanish and have opposite 
signs. As we shall see in the next section, it is these two last terms which are 
responsible for the particle’s drift. 


5 Driving force 


Substitution of Eq. (H51) for the collision times into Eqs. (1261) and (EZD yields 
the following expression for the average forces on two sides of the particle: 

{Fait)) = ± Eo I 1 -f Cl (uJaTa) tta ^a^ (t) COS Wat 

+ C2 (Wafa)^ tta SinUat + C 3 {iOaTaf a\ COS^ Wat |. (45) 

Here the sign of the right hand side is plus and minus respectively for the 
force on the left [a = 1 ) and right (a = r) sides of the particle, Fq is the 
magnitude of the force on each side for the passive regime given by OT, 

Fq = nksT = nmv^, (46) 


fa is the collision time for the passive regime for a molecule with the thermal 
speed vt = ^/ksT/m, 


= 'Taivr) 


2mvT 


and numerical coefficients are 



(47) 


(48) 
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The above expression (gSl) is a perturbation expansion of the average force 
up to second order in w, or more precisely, in the dimensionless parameter 

ea=UJaTa- (49) 

The approximation of zeroth order {Fa{t)) ~ ± Tq coincides with the force for 
the passive regime, Eq. m, and corresponds to the adiabatic approximation 
for the collision time which neglects the variation of the potential during 
the collision. As expected, in the adiabatic approximation, asymmetry of 
the static force magnitudes fa and of parameters of dynamic modulation 
(ocWq) do not show up, and the the total average force on the particle 
(F) = (Fi) + (Fr) is zero. 

The lowest order correction to the adiabatic approximation in linear in 
Wq (or Cq) and is presented by the second term in the right hand side of 
Eg. (1451) . If one considers a; as a variable responsible for perturbation from 
equilibrium, the truncation 

(F„(t)) « ±Fo I 1 + Ci (a;afa)aaCa^(t) COSWatj (50) 

can be interpreted as a linear response approximation. In contrast to the 
adiabatic approximation, asymmetry of static and dynamic parameters is 
manifestly present here, and the time-dependent part of the total average 
force {F) = (Fi) + (Fr) does not vanish. However, if in addition to the 
ensemble average one also evaluates the time average over the period of 
modulation, the function fa'^it) cosujat vanishes, and so do the total average 
force and the average velocity of the particle. The absence of the drift in this 
case can also be directly demonstrated solving the Langevin equation m- 
Thus the linear response approximation dlol) is insufficient to account for the 
operation of the system as a motor. 

The last two terms in the right hand side of Eq. (gSl) are of second order 
in oja and present the correction to the linear respond approximation. Both 
terms do not vanish under time averaging over the modulation period and 
thus can generate a drift of the particle. As one can check, the time averages 
of these two terms are of opposite signs. 

We finish this section with a remark that our derivation of Eq. (1451) for 
the average forces does not assume that the modulation amplitudes Oq, are 
small. An additional assumption Oa <C 1 and the linear approximation for 
the functions ^^"(t) 


fa^ ^ ^ ~ ^an sin (jJat ( 51 ) 

bring Eq. (|^5|) to order to the following form 

{Fa{t)) = ± Ed I 1 -I- Cl (uiaTa) CLa COSUJat — Cl (uJaTa) sin 2 Wat 

+ C2 (WaTa)^ Oa sinWat F {uJaTaY [^3 COS^ Wat — 3 C 2 sin^ Wq!] |. 

Clearly, only the last term, quadratic in both oja and Oa, contributes to 
the drift of the motor. Thus for the present model the drift is a nonlinear 
phenomenon with respect to both modulation parameters. 
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When both Wq and Oq, are small the drift is small too and numerical 
simulation becomes very time-consuming. In remaining sections we shall work 
with the more general expression (HSl) which does not assume the smallness 
of Oq,. 


6 Results 

In this section we discuss solutions of the Langevin equation with the 
average fluctuating force given by Eq. (HSl) . Although the equation is lin¬ 
ear, its explicit analytic solutions are rather bulky and not instructive. We 
therefore present solutions in a graphical form only. 

Let us first consider the case already discussed in Section II when the 
asymmetry is due to unequal magnitudes of the passive forces fi ^ fr, while 
the modulation frequencies and amplitudes for the left and right sides are 
the same, 


uJi= 0 Jr= 0 j, ai = ar = a, ^i{t) = ^r{t) = ^{t). (52) 

Introducing the parameter of asymmetry 6 by relations 

fr = Sfi, or fr = S~^Ti, (53) 

one obtains from 65]) the following expression for the total average force 

(F) = (Fi) + (Fr): 

{F{t)) = Fq |ci (wT/) (l — cosujt 

+ Cl (l — (5“^) sinwt 

-I- C 3 (wT;)^ (l — d~^) ^”"^(0 cos^ ujt |. (54) 

This expression is to be substituted into the Langevin equation (0, which 
we integrate numerically using the following time, space, and velocity units 

, 1 . mvT , xo 

to = -;^Ti = —^, xo = VTto, vo = — =vt- (55) 

2 Ji to 

The same units of course are employed in simulation. (Note that the accel¬ 
eration unit is oo = vo/to = fi/m, which means that molecules in the left 
interaction zone have acceleration of magnitude one). 

In dimensionless form the Langevin equation (HU) with the force (15411 and 
7 given by (7) reads as follows 

{V) = — 3 Cl h {V) + 2ci nauj {1 — S~^) ^~^(t) coswt 
dt 

-I-4 c2 (1 — d“^) ^“^(t) sinuii 

+ 4 C 3 n u}'^ {I — 6 ~‘^) (t) cos^wt, (56) 


where A = \/m/M, and a superposed tilde is used to denote dimensionless 
velocity V = V/vq, time t = t/to, frequency 6j = uj to, and the concentration 
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time 

Fig. 7 Average dis placement (x{t)) and velocity {v{t)) of the motor as functions of 
time for regime (E3 with parameters u) — 0.1, a — 0.5, and fr = fi = 1- Simulation 
data (solid lines) are averaged over about 2 • 10® trajectories. Theoretical curves 
(dashed lines) are solutions of the Langevin equation (1591) . 


of bath molecules n = nx^. Since the operations of taking a time derivative 
d/dt and an ensemble average (...) clearly commute, the equation for the 
dimensionless displacement X = Xjx^ follows from (1561) simply by the re¬ 
placement (F) = The solutions for {X{i)) and {V{i)) are presented in 

Fig. 2 for zero initial conditions by dashed lines. They are in good agreement 
with numerical simulation (solid lines) provided the modulation amplitude 
is not too high, a < 0.5. 

Motors with other types of asymmetry can be considered in similar ways. 
As another example, let us compare predictions of the theory and simulation 
results for the case when the left side of the particle is active, while the right 
side is passive, 

LOi = to 0, ai = a ^ 0 = Or = 0. (57) 

For this case, the general expression (HHll gives the following result for the 
ensemble-average fluctuating force on the particle 

{F{t)) = Fq |ci (wf) cosujt 

+ C 2 a (t) sinujt + C 3 (ujt)^ cos^wt|, (58) 

where all quantities are for the particle’s left side. The corresponding Langevin 
equation for the average velocity in dimensionless notations has the form 

(F) = — 3 Cl h (y) -|- 2 Cl h a w (t) cos Cot 

dt 

+ Ac 2 naui'^ sin Cot + 4 C 3 X"^ fia^ Co^ cos^wf. (59) 

The solution of this equation and of a similar equation for the dimen¬ 
sionless displacement {X{t)) are presented by dashed lines in Fig. 7. Again, 
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within the range of its validity, lot 1, the theory agrees with simulation 
provided a < 0.5. As the variation amplitude a is further increased, the the¬ 
ory progressively overestimates the drift velocity. We believe this is due to 
neglected effects of the variation on the dissipative force in the Langevin 
equation. A small offset between theoretical and experimental curves visible 
in both Fig. 2 and Fig. 7 can probably be attributed to non-Markovian effects 
which are also neglected in the presented theory. We hope to address these 
issues in future work. 


7 Conclusion 

The relations between internal conformational dynamics of complex biomolec- 
ular systems and their ability to move directionally are of considerable in¬ 
terest in many fields. Ab initio simulation of coupled conformational and 
translational dynamics is usually of high computational cost and often not 
illuminating. In this paper we studied a reduced model where the internal 
dynamics of a motor is not considered explicitly but is assumed to result 
in periodic variation of the strength of the microscopic potential through 
which the motor interacts with molecules of the surrounding thermal bath. 
The model is simple enough to allow the evaluation of the driving force in 
an analytic form. In contrast to many models with variation of external pa¬ 
rameters, our model shows no directional motion in adiabatic approximation 
of infinitely slow (reversible) variations. Moreover, the leading-order correc¬ 
tion to the adiabatic approximation (linear in the variation frequency to) is 
also insufficient: Without the last two terms quadratic in to, the Langevin 
equation (IHIH) . or (1^ . has a no-drift solution. Interestingly, those two terms, 
if only one of them is left in (j59l) . generate the drift in opposite directions. 
Thus, being the interplay of two terms, the direction of motion of the motor 
can hardly be predicted with qualitative arguments. 

The model can be extended in many ways. Since the role of internal vari¬ 
ations is to drive the motor out of equilibrium, it is clear that a strictly 
periodic and deterministic character of variations is unnecessary to main¬ 
tain the drift. One may expect a similar mechanism of directional transport, 
for instance, for motors driven by internal conformational transitions which, 
due to the fuel consumption, do not satisfy the equilibrium detailed balance 
condition. Models with stochastic transitions between two and more inter¬ 
nal conformational states were discussed earlier [21], but they implied the 
existence of the external potential(s), whereas our motor is autonomous. 


Acknowledgements I thank to G. Buck, G. Parodi, and J. Schnick for discus¬ 
sions, and to an anonymous referee for valuable comments. 


Appendix 

This Appendix (not included in the published version) presents perturbation 
solutions of Eq. (41) to second order in the small parameter e = lot^ . The 



19 


equation has a form 


1 + cix + ec2X^ + e'^ csx^ = 0 


where x = t* jr^ and 



Substituting into the equation the second-order ansatz 


X = xo + e xi + X2 


and discarding terms of order higher than two yields 

0\ /O 

1-I-Cl(xo-b e xi-I-e X 2 )-b ec 2 (Xq- b 2 exo Xi)-b e C 3 Xq = 0. 
Equating to zero contributions of each order separately 


1 -b Cl Xo = 0, Cl X1 -b C2 Xq = 0, C1X2 -b 2 C2 Xq X1 -b C3 Xq = 0, 


one finds 



1 


Xo = - 


Cl 


For T* = XT° these relations lead to the solution dHI). One can show that 
two other solutions are singular (diverge as e —> 0) and have no physical 
meaning. 
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